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I.  INTRODUCTION 


Virtual  cathode  formation  in  vacuum-propagating  beams  may  be 


more  widespread  than  previously  recognized.  The  '"spin  death-"  of 


parallel  beam  energy  by  ErxB?  rotation  in  front  of  a density  pileup. 


may  assit  or  steepen  the  axial  density  gradient,  since  the  rotation 


speed  tends  toward  Brillouin  flow,  u)=u  /2.  It  may  be  that  the 

c 


2 2 

propagation  limit  (on  the  beam  density)  is  not  u>  'Vj  in  the  applied 

P c 

2 

or  total  magnetic  field,  but  rather  w (in  the  presence  of  self- 

p max 


field  pileup)  small  enough  that  cE^/B^Vq  (Vg=beam  electron  speed), 


with  given  self-consistent ly  from  the  charge  distribution. 


In  what  follows,  self-consistent  BGK-like  axially  dependent 


equilibria  will  be  described,  in  which  the  space  charge  effects 


determine  the  axial  extent  and  profile  of  the  beam.  Generalization 


to  axially  nonuniform  magnetic  fields  is  relatively  straightforward. 


Sections  III  and  IV  treat  the  problem  in  a purely  one-dimensional 


way,  keeping  the  electrostatic  effects  of  E electric  fields  but  not 


the  effects  of  beam  rotation  caused  by  Er.  Section  V generalizes 


these  results  to  include  E . Section  VI  summarizes  the  conclusions 

r 


of  this  introductory  work. 


II.  APPROXIMATIONS  AND  ASSUMPTIONS 


To  calculate  the  electron  density  n(s)  at  a point  s along 


magnetic  field  lines  in  terms  of  an  effective  potential  U(s)  for  the 


parallel  motion,  it  is  assumed  that  all  changes  of  U(s)  and  n(s)  with 


time  occur  much  slower  than  the  time  required  for  a typical  particle 


transit  from  s=0  to  the  end  of  the  system  and  back,  and  that  the  beam 


is  on,  supplying  electrons  at  a constant  rate,  for  a time  much  longer  than 


the  factor 


✓ 2 9 2 2 

Y = T1  + (pl(  + pP/iDqC  , 


which  occurs  in  the  n(s)  integrals. 

The  collisionless,  magnetic-moment-preserving  behavior  assumed  here 
implies  that  instabilities  do  not  disrupt  the  electron  motion.  This 


t 


requirement  may  be  satisfied  since  the  virtual  cathode  formation  will 
cause  inhomogeneity  of  the  density  on  the  scale  of  one  beam  Debye  length. 
Thus,  even  though  reflection  of  a beam  with  small  parallel  momentum 
spread  generally  allows  counterst reaming  instability  in  homogeneous 
plasmas,  this  instability  may  not  be  able  to  grow  to  significant 
amplitude  in  the  space  provided,  all  the  more  because  density 
gradients  absorb  the  unstable  waves.  Radial  motion  due  to  E^xB^ 
is  neglected.  This  additional  self-pinching  due  to  E must  be  offset 
by  centrifugal  effects  arising  from  the  curvature  of  the  helical  magnetic 
field  lines. 


III.  DENSITY  IN  TERMS  OF  POTENTIAL 


The  density  at  s is  comprised  of  forward-going  and  reflected 
electrons : 


n(s)  = n+(s)+n_(s) 


The  velocity  distribution  of  the  forward-going  portion  is  specified 
at  s=0: 


f+(vl  , vA,s»0)=N<5(  | v|-v())H(v)H(  Iyq-VjJ) 


where  H(x)  is  the  Heaviside  function:  H(x)=l  for  x>0,  =0  for  x<0,  and 


I 


y^  Is  a maximum  value  of  v^.  This  distribution  is  nonzero  on  an  arc 
of  a circle  in  v space  at  s=0.  See  Fig.  1. 

The  velocity  distribution  of  the  reflected  electrons,  f_(v|(  ,vx,0) 
follows  from  reflection  conditions  and  so  is  not  specified  a priori. 


One  can  now  use  the  integrated  moment  form  of  the  one-dimensional 


steady  state  Vlasov  equation  to  express  n(s)  in  terms  of  f+(v,0)  and 
4>(s)  : 


n(s)  = 2h[B(s)/B(0) ] 
x = v , y = v,  , 


ydy 


xdx(x^-x^)  ^^f(x,y,0) 

r s 


(3) 


where 


2 

>0 


mr  1 


(m 


is  the  square  of  the  'escape'  velocity  required  to  reach  point  s, 
and 

l = Sup  x , (y) 


(5) 


s <s 


is  the  maximum  value  taken  on  by  x , at  points  between  s'=0  and 


s =s . 


Particles  with  midplane  parallel  velocity  x>xg  contribute 


to  the  density  at  s;  those  with  x<x^  or  ) have  been  reflected 
at  some  s' <s  and  do  not  appear  at  s.  Because  of  the  beamlike  nature 


of  f+(v,0),  there  are  no  electrons  with  x less  than  some  Xq  (i.e.. 


with  y greater  than  some  y^) ; thus  for  s less  than  some  Sq  there  are 


no  turning  points.  For  s>Sq  there  are  turning  points  until  s >M.  The 


point  sQ  is  specified  by  its  potential,  i.e.. 


with  x given  by  Eq.  (4)  with  s replaced  by  s_.  The  minimum  parallel 
s0  u 


velocity  at  the  reference  plane,  xQ,  is  a parameter  of  the  distributi 


funct ion. 


For  points  s>sq>  the  lower  boundary  of  the  x-integrat ion  in 
Eq.  (3)  is  E,  as  given  in  Eq.  (4).  For  s<sq»  the  minimum  is  simply 
Xq.  In  both  cases  the  maximum  can  be  taken  as  Vq,  as  is  clear  from 
Fig.  1.  Thus,  one  has  when  B(s)rB(0) 


fys  , rV  0 

n+(s)  = y N . d (y  ) 


J0  ‘2  2/2  2 
7.  or  x_  /x  -x 
0 s 

= 2ttV  n/v'T-x*'  for  s>s  , 

0 0s  0 


d (x  ) 2 

°(*x  +y  -vQ) 


= 2nV0N  [/v^-x;  - /v^-y^]  for 


s<s0  ‘ 


In  Fig.  2,  the  forward-going  density  n+(s)  is  shown  vs  xg,  i.e.,  vs. 
the  potential  difference 


A <f  (2e/m)  [*(0)-»<s)  ] 


The  density  of  reflected  particles  is  calculated  similarly. 

One  assumes  there  is  an  s=M  such  that  no  reflections  occur  beyond  M. 

Then  n (s)  is  comprised  of  particles  that  go  past  s but  not  past  M: 

...  2 2 2 


n_(s)  = y N 


Mln  V>’0  . rX m , . 2. 

d (y  ) X 6(/x+y  -vn) 

2 >2/22 

y X tx  -x 


Note  that 


/ 2 2 

2iV  NvV-x  for  s <s<M 
0ms  0 


. „ / 2 2 / 2 2,  r 
= 2ttV„N[/x  -x  -»x.-x  1 for  s<s_ 
0 m s 0 s 0 


x"-x2  = (2e/m)  ( <t> ( M)  — <J> ( s)  1 
m s 


The  density  at  s is  then 


(14) 


I j 


L 


n(s)  = n+(s)+n_(s)  = 2TrNvfl  [ /v2-Af  - 


+ /A/  -Af  - A2-^} 
m U 


with  any  radical  dropped  once  it  becomes  imaginary.  For  convenience. 


2 

let  i^A^P/v  The  physical  range  of  the  "independent  variable"  1J1 


is  then  from  zero  (i p(0)-0)  to  ip  (<P(M)S<P  )•  The  density  n(s)  is 

M m 

shown  vs.  ip  in  Fig.  3. 


IV.  SELF  CONSISTENT  SOLUTION  FOR  AB=0 


The  equation  for  n(s)  in  terms  of  ip ( s ) is  now  coupled  to  a 
Poisson-like  equation  for  i|'(s)  whose  right-hand  side  is  n(s). 
The  Esulting  equation  has  the  form 


a^U1)  = w^(il,3^/ar) 


(15) 


where  u is  a functional  of  tp  (and  in  general,  its  derivatives) 
P 


and^£is  the  Laplacian  operator.  The  simplest  model  is  a purely 


one-dimensional  one  with  |J’(s)=<)> (s)  (i.e.,  purely  electrostatic  potential) 
« 2 2 

and^£=d  /ds  . In  terms  of  the  scaled  distance  z ^s • ( 2co^Q/Vy)  the 
equation  is 


d2<p 

dz2 


[ 2 \p+  /iji  -\p  ] 

U m 


n-2^r+^  i 

u m 


(16) 


2 2 

where  xq^vq  ant^  where  the  normalization  N has  been  evaluated  in  terms  of 


2 2 

the  density  at  s=0  and  subsumed  into  the  definition  of  z via  u>  Q-^llnQe  /m. 


It  is  not  known  whether  this  nonlinear  equation  always  has  unique 
solutions.  When  the  spread  in  parallel  velocities  of  the  beam  is  small 


we  can  cast  this  equation  in  a soluble  form  if  all  electrons  are  reflected 


1 2 

somewhere:  then  =1,  i.e.,  e I ♦(O) -'J’(M)  ] = — mv_,  and  the  right-hand  side  of 
m Z U 


(16)  is 


For  small  angle  spread  we  may  write  i(/0=l-c  (see  Fig.  1)  where 

2 

in  terms  of  the  pitch  angle  spread  0^,  e=6q/2.  Expanding  in  e, 
the  right-hand  side  becomes 


a 

/l~? 

where  a->l  as  e-K). 

The  resulting  equation. 


dz^  /l-ii i 


is  a turning  point  problem  and  can  be  solved  by  quadrature: 


z = 


_1_ 

3a 


[(I-/I-4.)  (1+2 ✓!-</) 


1/2 


1 


(17) 


(18) 


giving  ^(z)  implicitly  as  z(i}j).  In  particular,  the  condition  that 

no  particles  are  reflected  beyond  some  point  s=M  (and  that  all  are 

reflected  somewhere)  gives  t)'  or  M,  subject  to  <1: 

m m 

m - -f-  Hi-/iT)(i+2/iT)1/2l  • U9) 

3a  m m 


See  Fig.  4. 

In  the  limit  of  no  angle  spread,  all  the  electrons  are  reflected 
at  this  same  point.  One  can  then  use  Eqs.(14)  and  (18)  to  obtain  n(s) : 


n(s)  = l/u(s)  u(s)  = /1-iJj  (20) 

u3-]u2-u+|  (1-9Z2)  = 0 (21) 

with  z»s*2w  „/v _ as  before. 
pO  0 

In  the  case  that  not  all  electrons  are  reflected  (i.e.,  U(s)  has 
2 

a maximum  less  than  mv^/2)  the  one-dimensional  electrostatic  problem 


can  be  similarly  reduced  to  numerical  quadrature: 


,2 , __ 

dz“ 


with 


3=[l-2/^+>>ni]_1  . 


The  solution  is  obtained  by  numerically  integrating 


(22) 


(23) 


dz  = 


14' 


where  ^ 

W(i|<)  = a 


21W  -W(*)] 


1/2 


di^'l  </l-ijwV\jj  — 2 U *1 

m U 


3/2  „ , 3/2 


,3/2 


3/2. 


3/2, 


(24) 


= ll+v  -2^  — ( 1 — ti>)  — ( tl'  -*)J'N-2(*  -*)-"‘]  (25) 

1 m o m U 

2 2.  2 

Here,  as  before,  ^o=x0^vO  t^e  COS  t^le  *3eam  angle  spread;  and 

^max  t^le  max3rauin  value  of  W in  the  range  of  integration  over  iji. 


V.  EFFECT  OF  RADIAL  ELECTROSTATIC  FIELD 


Although  the  treatment  in  Eqs . ( 16 ) —(25)  is  strictly  one-dimensional, 
i.e.,  involves  only  electric  fields  along  s,  a formal  generalization 
is  also  possible  when  electron  parallel  streaming  energy  is  converted 

— V -f 

to  azimuthal  Eq*B  rotation  due  to  the  radial  electrostatic  self-field, 
EqCt.s).  Since  any  increase  in  the  energy  of  this  azimuthal  motion, 

\ m(cEQ/B)2  , 


must  come  from  the  initial  streaming  energy,  it  can  be  treated  as 
part  of  the  effective  potential  for  parallel  motion,  U(s) , which  now 
becomes 

U(r,s)  = ™ y2B^B(s)-e4!(s)  +|  C2  [EQ(r  , s) /B(  s)  ]2  . (2 

As  mentioned  in  Sec.  I,  this  can  be  written 

U(r,s)  = y2B”1B(s)-ei>(r ,s)  (2 


where  4>  now  includes  all  effects  due  to  the  electric  field,  including 
the  radial  field,  and  where  the  magnetic  moment, 

m 2-1 
^2y0B0  * 


is  assumed  conserved. 

For  arbitrarily  large  beam  density,  azimuthal  current  now  makes 
B(s)  nonuniform: 


and  the  natural  coordinates  ( M j.  ?)  become  noncylindrical , complicating 
2 

the  V operator.  However,  to  demonstrate  the  physics  of  the  electrostatic 
effects  alone,  we  may  now  restrict  the  remaining  discussion  in  this 
section  to  the  case  of  small  heam  rudius,  low  density,  or  strong 
Writing  Eq . (28)  in  dimensional  form  with  S(r ,s)*?^+^s(r ,s) , one  can 
show  that 


V me  / 

' % 2 

pO  \ me  / 

with  the  beam  radius  and  the  plasma  frequency  based  on  the  beam 

density  at  s=0.  Thus  in  Eqs.  (30)-(35),  we  shall  require 

2 

r,  u> 

— « 1 ( 3C 

C <0 

c 

so  that  the  magnetic  self-fields  are  weak,  while  the  electrostatic 
self-fields  are  still  strong  enough  to  be  important.  In  this  case, 
because  of  the  electrostatic  nature  of  the  problem,  one  has 


EQ(r,s)  = - --  <p  (r  ,s) 


and  straight  magnetic  field  lines.  Thus  we  can  write 


>(r,s)=<j>(r,s)-(mc2/2eB“)(^-) 


for  the  case  B(s)=const.=Bp. 

Because  reflection  of  electrons  depends  only  on  the  potential  for 
parallel  motion,  an  azimuthally  symmetric  nonneutral  distribution  of 
electrons  will  have  its  density,  n(r,s)  given  by  Eq.  (3),  as  before, 
except  for  the  modification  in  the  definition  of  $ . For  the 
distribution  function  of  Eq.  (2),  the  results  of  Eqs.  (6)-(14) 
remain  valid.  A self-consistent  equilibrium  for  straight  magnetic  field 
lines  can  now  be  obtained  in  the  form  of  Eq.  (15),  letting 

^-^yMr.s)  (33) 

mv0 

be  the  scaled  electrostatic  potential,  and  defining 

[v0/(2eB0/mc) ]2~b  , (34) 

2 

we  may  simply  replace  ip  by  'f-t-b ( / 3 r ) in  the  right-hand  side  of  Eq  . 

2 

(16).  Note  that  b is  just  1/4  of  the  (gyroradius)  , so  that  the 
natural  radius  scaling  is  of  the  order  of  the  gyroradius.  Eq . (16) 
then  becomes,  for  the  2-dimensional  electrostatic  problem. 


7 t?  (r  If)  + ^ = -{/i-*-bO*/ar) 

az  

-2/cos20„-il)-b(a 


0-i|)-bO<jj/ar)  ‘ 


h/lax^  (tj,+b(aij</ar)2  ]-[il/+b(ail,/ar)2]  ' 
t |l-2coseQ+^ilaxz  [ili+b(alj)/ar)^T | 


where  0.  is  the  initial  spread  in  pitch  angles  of  beam  electrons,  and 
where  the  operation  Max^  is  taken  at  fixed  r.  The  nonlinear  partial 
differential  equation  (35)  appears  suitable  for  moderately  simple 
computer  solution  by  overrelaxation  methods. 

The  condition  (30),  together  with  the  requirement  that  the  electron 
gvroradius  be  small  compared  with  the  beam  size  r^,  give  the  requirement 

vo/c<<<“c/up)2  * 

where  to  eB../mc. 
c 0 


VI.  CONCLUSIONS 


L 


In  the  one-dimensional  case  ol  purely  electrostatic  repulsion 

(£^>0,  11^=0) , self  reflection  of  the  entire  electron  population  occurs 

in  a finite  distance  even  if  there  is  no  pitch  angle  spread.  For 

small  or  zero  spread,  this  virtual  cathode  occurs  at  a distance  L^v„/6ui  „ 

'v  0 pO 

i.e.,  at  less  than  one  beam  Debye  length  based  on  the  initial  kinetic 
energy  per  particle.  For  a cold  beam,  all  the  beam  electrons  are 
reflected  at  the  same  point,  namely  ;it  this  distance  L.  The  density 
profile  is  then  given  by  n(s)=l/u  with  u the  solution  of  the  cubic 
equation  (21);  this  profile  is  shown  in  Fig.  5.  The  purely  electrostatic 
one-dimensional  problem  with  angle  spread  has  been  reduced  to  quadrature, 
as  given  by  Eqs.  (24)  and  (25). 

For  the  case  with  F.^0,  the  scaled  electrostatic  potential 
^=2e$/ mv^  is  given  by  the  parabolic  equation  (35).  This  equation  is 
suitable  for  moderately  simple  computer  solution  by  the  method  of 
over re  1 axat ion . 


FIGURE  CAPTIONS 


RffenCMC** 

Fig.  1 kU^plane  velocity  space,  x'v)(  , y*v1,  at  spatial  point  s = 0. 

Beam  input  velocity  distribution  is  nonzero  on  the  spherical 

cap  with  x>x_..  Beam  electrons  with  x_‘-x<x  are  reflected 
U Us 

somewhere  in  the  interval  (0,s)  and  do  not  contribute  to  n(s). 
Reflected  particles  (x<0.  also  contribute  to  n(s).  There  may 
be  an  x^<Vq  such  that  electrons  with  x^<x^v(|  are  never  reflected, 
if  In  this  case,  -v  <x<-xvl  does  not  contribute  to  n(s). 


Fig.  2 


Fig.  3 


Forward-going  electron  density  n+(s)  at  s [Eqs.  (7)  and  (8)], 
vs.  the  potential  difference  [Fq.  (9)]  between  point  s and 
s»0;  shown  for  beam  pitch  angle  spread  6^  at  s=0  given  by 
sln260;y2/v2-0,  0.1,  0.3,  0.5. 

Total  electron  density  n(s)  at  s,  vs.  the  normalized  potential 
difference  between  point  s and  s*0  [Eq.  (14)].  Shown  for 
iJ'q'O.?,  i^»0.95.  The  lower,  dashed  curve  describes  the 
never-reflected  electrons  at  s>M,  i.e.,  beyond  the  potential 

maximum. 


Fig.  4 ii»(z)  for  z<M,  from  Eq.  (18). 
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